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ABSTRACT 

We investigate the relationship between the star formation rate per unit area and the 
surface density of the ISM (the local Kennicutt-Schmitt law) using a simplified model 
of the ISM and a simple estimate of the star formation rate based on the mass of 
gas in bound clumps, the local dynamical timescales of the clumps, and an efficiency 
parameter of around e « 5 per cent. Despite the simplicity of the approach, we are 
able to reproduce the observed linear relation between star formation rate and surface 
density of dense (molecular) gas. We use a simple model for the dependence of II2 
fraction on total surface density to argue why neither total surface density nor the Hi 
surface density are good local indicators of star formation rate. We also investigate 
the dependence of the star formation rate on the depth of the spiral potential. Our 
model indicates that the mean star formation rate does not depend significantly on 
the strength of the spiral potential, but that a stronger spiral potential, for a given 
mean surface density, does result in more of the star formation occurring close to the 
spiral arms. This agrees with the observation that grand design galaxies do not appear 
to show a larger degree of star formation compared to their fiocculent counterparts. 
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1 INTRODUCTION 

1.1 Observational Background 

The dependence of the star formation rate per unit area on 
the surface density of the ISM in a galaxy is one of the most 
highly speculated problems in extra-galactic surveys, as well 
as theoretical analysis of star formation. The star formation 
law predicts both the amount of star formation and the de- 
gree of stellar feedback in a galaxy, and is consequently a 
vital requirement for models of galaxy evolution (e.g 



Tan, 



Silk & Balland |1999| |van den Bosch||2000| |Springel||2000| 


Hou, Prantzos & Boissier|2000 


Springel & Hernquist|2003 1. 



Schmidt ( 1959 1 related the star formation rate per unit 
volume, psFR, to the gas density, p, in a galaxy accord- 
ing to psFR tx p™. An observational law for star forma- 
tion rate per unit surface area, Esfr in terms of the mean 
gas galactic surface density, E, viz. Esfr oc E" was estab- 
lished by Kennicutt] |l989), and since then numerous sur- 
veys have attempted to find a universal value of n. |Kennicutt| 
( 1989 1998 1 investigated a global star formation law in star- 



forming galaxies and found n — 1.3 ± 0.3 from a sample of 
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97 galaxies including normal spirals and starbursts, for gas 
surface densities of a few M© pc~^ to around 10** Mq pc~'^, 
although it should be noted that at a given E, the scatter 
in log Esfr is generally more than ±0.5. Other surveys in- 
clude Gao & Solomon (2004b), who trace both the densest 



gas, using HON, and the total gas in a galaxy. They obtain 
different relations depending on what is meant by E. For E 
measuring dense gas they find n=l whereas for E measuring 
total gas density they find 1.73. 

The local relationship between star formation rate per 
unit area and surface density (what we call here the local 
Kennicutt-Schmidt law), has also been examined for indi- 
vidual galaxies. [Wong fc Blitz| ( |2002[ ) used radially averaged 
values of the star formation rate and surface densities for 6 
galaxies, finding n=l. 2-2.1. Kennicutt et. al. 2007 also ob- 
tained n = 1.56±0.04 for M51, by sampling the Ha emission 
and total surface density over 500 pc size regions (see also 



Schuster et al. 2007 1. There is a slightly less steep correlation 
(n = 1.37±0.03) for the surface density of molecular gas, but 
no correlation with the HI gas. Similarly Heyer et al. (20041 



found n = 1.36 for M33, when considering the molecular gas 
alone, but a much steeper dependence (n = 3.3 ±0.1) on the 
total gas surface density. 
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In fact the general consensus from these results is that 
a well-defined local star formation law holds only for the 
molecular gas in these galaxies. This conclusion is further 



endorsed by recent results from the THINGS survey ( Bigiel 
et al.|2008 1 , which show that there is no correlation between 
the Hi and the star formation rate - Espr is multivalued for 
a given density of Hi. Their results also indicate a fairly 
sharp transition from a regime where gas is predominantly 
Hi, to where gas is mainly molecular. Thus in the Hi regime, 
SsFR varies very steeply with E, whereas in the H2 regime, n 
is roughly 1. For the total gas surface density, n lies between 
1 and 3, although it is evident that the data does not give 
a good linear fit. 



1.2 Theoretical interpretations 

Several theoretical explanations of the Kennicutt-Schmidt 
law have also been advanced. All interpret E as being the 
total surface gas density. Since star formation occurs in a 
chaotic or turbulent environment, a natural explanation is 
that turbulence somehow regulates how much gas exceeds 
the high densities required for star formation. |Elmegreen| 
1 2002 1 and Krumholz & McKee ( 2005 1 assume a probability 



density function of densities in a turbulent regime to obtain 
a star formation law with n — 1.4. For jElmegreen (20021, 



this involves assuming an unknown star formation efficiency 



(e), but Krumholz & McKee (20051 instead determine the 



star formation rate efficiency, thus eliminating e. This is es- 
sentially the derived star formation rate divided by the total 
possible star formation rate if all gas at a particular density 
was turned into stars. 



Krumholz & Tan (20071 and Krumholz & Thompson 



(20071 further compare the turbulence regulating model 



of star formation with observations. IKrumholz fc McKeel 



( 2005 1 predict that the star formation rate efficiency is 



~ 0.013, independent of density. Their observed estimates 
of the star formation rate efficiency ( [Krumholz fc Tan|2007[ ) 
are consistent with this prediction. However the uncertain- 
ties (in particular for lifetimes for a given tracer) do not rule 
out a star formation rate efficiency which increases with den- 



sity either. Krumholz fc Thompson (20071 also determine 



CO and HCN luminosities, and show that the star forma- 
tion rate is linear with relation to L(HCN) and nonlinear for 
L(CO). 

An alternative possibility is to estimate the star forma- 
tion rate from cloud-cloud collisions. In this case Esfr is 
proportional to the surface density over the collision time, 
which in turn is dependent on the shear in the disc, mea- 
sured by the orbital angular velocity Q,. This yields an al- 
ternative form of the Schmidt law, Espr oc eQ,Tj^ where e is 
the star formation efficiency and N = 1 ( Wyse|[l986 Wyse 
|fc Silk|p89l [Tan|[2000| . [Silk| ( (T997| explicitly includes the 
star formation efficiency by assuming that star formation is 
self regulating in the disc, obtaining ^sfr oc E^'^®. 

In addition to the observational surveys, several nu- 
merical simulations have investigated star formation rates 
in galaxy simulations ( Li et al.|2006 Robertson & Kravtsov 



|2008[|Tasker fc Bryan|2008[ ). However a major problem with 
numerical evaluations of the star formation rate is that it is 
dependent on what assumptions are made about the con- 
ditions for forming sink particles at high densities, or for 
taking account of possible forms of stellar feedback. A den- 



sity threshold is required, as well as a star formation efh- 
ciency, and so the degree of star formation in these simu- 
lations is also dependent on both these parameters. Saitoh| 
et al. ( 2008 1 performed simulations with varying thresholds 



and star formation efficiencies, concluding that the star for- 
mation rates and Kennicutt-Schmidt law do not significantly 
change provided there is a high threshold density. However 
an estimation of the star formation rates from first princi- 
ples, without these approximations, is beyond current capa- 
bilities. 



Tasker & Tan ( 2008 1 perform similar models to those 



of the current paper, and likewise do not include any star 
formation prescription. In their simulations there is no un- 
derlying spiral potential, rather numerous flocculent spi- 



ral arms cover the disc. In agreement with Dobbs ( 2008 \ 
they find that GMCs form by agglomeration via collisions, 
although collisions occur throughout the disc whereas in 



Dobbs (20081, collisions between clouds are largely confined 



to the spiral arms. They do not make any estimates of the 
star formation rate from their results, however they do show 
that the number of collisions between clouds is consistent 
with the model presented in |Tan| (j2000| . 

1.3 The current paper 

One approach to trying to understand the physics under- 
lying the basis of the Schmitt-Kennicutt relations is to un- 
dertake ever more detailed numerical simulations with ever 



increasing quantities of input physics (e.g. Susa 2008 Shetty 
fc Ost rikcr 2008| |Agertz et al.|2009| . In this paper, we take 
an alternative approach and ask the question: how little in- 
put physics, and how few assumptions do we have to make, 
in order to obtain relations which resemble the observational 
findings to a reasonable degree? In this manner we hope to 
be able to obtain some understanding of what the funda- 
mental drivers for such a relationship might be. 

The calculations presented in this paper make use of the 
simplified numerical simulations which have already been 



published in Dobbs ( 2008 1 . The simulations were performed 



using Smoothed Particle Hydrodynamics, and include mag- 
netic fields and self-gravity. We do not include the star for- 
mation process itself nor any subsequent stellar feedback. 
Thus the gas is supported against collapse by thermal and 
magnetic pressure in the lower surface density calculations. 
In the higher surface density results, we run the calculations 
until collapse occurs. The difficulty with the assumed lack 
of feedback is that we are restricted to relatively low surface 
density calculations. The calculations which have been se- 
lected for analysis in this paper are shown in Table [T] These 
calculations, with the exception of model C, are described 



in Dobbs (20081, but we also provide details below. 



1.3.1 Details of numerical models 

The calculations model a 3D gaseous disc between radii of 
5 and 10 kpc. The gas is assumed to orbit in a fixed galactic 



gravitational potential. The potential includes a halo ( Cald- 
[well fc Ostriker|[T98T| , disc ([Binney fc Tremaine]|1987| ) and 



4 armed spiral component ( Cox fc G6mez|2002 K The gas is 
initially assigned velocities according to a flat rotation curve 
determined by the disc potential, and in addition a velocity 
of dispersion of 6kms~'^ is superposed. 
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Model 


S (Mq pc-2) 


F (%) 


Qc 


Qh 


A 


4 


4 


0.5 


5 


B 


8 


4 


0.25 


2.5 


C 


16 


4 


0.25 


2.5 


D 


20 


4 


0.1 


1 


E 


4 


2 


0.5 


5 


F 


4 


8 


0.5 


5 


G 


4 


16 


0.5 


5 



Table 1. The calculations from pobbsl | |2008[ l which are used in 
this paper. All the calculations use 4 million particles and are 
isothermal. In these calculations, half the gas cold (100 K) and 
half warm (10* K), and the ratio of thermal to magnetic energy 
(/3) is 0.4 in the cold gas. F is a measure of the strength of the 
potential, and Q is the Toomre parameter (Qc for the cool gas and 
for the warm gas). The determination of all these quantities is 
described in |Dobbs| (2008^ . Models A-D compare surface density, 
and E-G shock strength. 



Models A-D have identical potentials but vary in their 
mean gas surface densities and therefore compare surface 
density. The total mass of the disc is 1 x l(f Mq in the 
4M0 pc"^ model, and 5x lO'' Mq in the 20 Mq pc"^ model. 
The simulations all use 4 million particles, hence the highest 
mass resolution is 250 Mq, and the lowest 1250 Mq. Models 
A, E, F and G have the same mean gas surface densities but 
adopt different strength for the spiral potential. 

In keeping with our aim for minimal input physics, the 
simulations are very simplistic. In particular, all the cal- 
culations in Table [T] assume an interstellar medium which 
has two isothermal components, one cool and one warm. 
We omit thermal considerations and so there is no transi- 
tion between the two phases; the cool gas remains cool and 
the warm gas remains warm, throughout. We use the same 
thermal distribution in all the calculations and only vary 
the global surface density and/or the shock strength. The 
cool gas is taken to have a temperature of T = 100 K, and 
we will think of it as representing molecular gas (H2). The 
warm gas is taken to have T = 10* K, and we will think of it 
representing atomic gas (Hi) . In all cases the cool and warm 
gas comprise equal mass in the simulations. It is of course 
possible to include thermal effects such as heating and cool- 
ing - see, for example ( jPobbs et al.||2008J (although these 
did not include self-gravity) but that is not the purpose of 
the current exercise. 

The initial scale heights of the warm and cold warm 
components are 150 and 400 pc respectively, giving a mean 
smoothing length of 40 pc. However with time the scale 
heights decrease to 20-100 pc and 300 pc. 

These calculations do not include sink particles, al- 
though the gas is self-gravitating. In the lower surface den- 
sity results (E 5C IOMqpc"'^), the cool gas is sufficiently 
supported by magnetic and thermal pressure, and the (typ- 
ically supersonic) velocity dispersion of the gas. However 
in the higher surface density results, runaway gravitational 
collapse does occur. The calculations were run for either 
250 Myr, or until the calculation is halted by collapse. The 
maximum timestep in the calculations (and frequency of 
dumps) is 2 Myr, but individual particles timesteps can be 
much less ( |Bate||1995 1. 



2 ESTIMATING Esfr IN THE NUMERICAL 
SIMULATIONS 

In this paper we shall assume that star formation occurs 
only in those regions of the ISM which are bound, and at a 
rate determined simply by the local dynamical timescale. 



2.1 Gravitationally bound gas 

We must first find a means of identifying which regions of 
the ISM are gravitationally bound at a particular moment. 
It needs to be borne in mind that while unbound gas can 
obviously become bound, it is also possible for bound gas 
to become unbound, despite the lack of feedback. This can 
come about for example as the gas accelerates out of a spiral 
arm, flowing over the ridge of the spiral potential, becoming 
sheared and longitudinally stretched as it does so. 

We determine the mass of bound gas using the output 
from a simulation at a given time frame (after 250 Myr ex- 
cept for models C and D, where the times are 200 Myr and 
140 Myr). The particular timeframe selected is not impor- 
tant, provided the time is not too near the beginning of the 
simulation, where the amount of bound gas will be underes- 
timated. 

In order to determine the mass of bound gas, we first 
sort the particles according to density. We select the most 
dense particle, and all particles within a radius of a smooth- 
ing length of that particle. Then we determine the gravita- 
tional, kinetic, thermal and magnetic energies for this group 
of particles. If {Ek + 2Eth + Emag)/Egrav < 1, the gas is 
assumed to be bound and the radius increased until the gas 
becomes unbound. Then the particles which are bound are 
recorded and their mass added to the total mass of bound 
gas. If on the other hand the gas is unbound, the gas par- 
ticles are discarded from the list. We also performed this 
analysis using solely either just the kinetic energy or both 
the kinetic and thermal energy. 

Using this method, the bound gas is situated in discrete 
clumps (which are spherical by assumption). Fig[l]shows the 
location of bound gas for a small section of the disc, for the 
simulation Model B with E — SMgpc"^. The colours indi- 
cate whether just the kinetic, the kinetic and thermal, or the 
kinetic, thermal and magnetic energies are included. As ex- 
pected, when all three energies are included, the extent of a 
clump which is bound becomes smaller, whereas the clumps 
are more extended (particularly true in the higher surface 
density calculations) when only the kinetic energy is consid- 
ered. In other simulations which investigate the Kennicutt- 



Schmidt law, magnetic fields are not included ( 


Bottema 


20031 ITasker & Bryan||2006| 


Li, Mac Low & Klessen||2006| 


Booth, Theuns & Okamoto 


20071, so by default only the 



kinetic and thermal energies can be used. 

In Fig. [1] the bound gas is dominated by two massive 
clumps which are just leaving the arms. These are more ex- 
tended since they have higher densities and lower velocity 
dispersions than clumps in the spiral arms. The gravitational 
energy of these clumps is also high since they are centrally 
condensed. Essentially these might be taken to represent 
massive Giant Molecular Clouds (GMCs), which formed in 
the spiral arms and are now entering the interarm regions, 
leaving only less massive clumps in the spiral arms. We also 
see that there are numerous bound clumps in the interarm 
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Figure 1. The location of bound gas is sliown for a small section 
of the disc from the E = 8M0pc~^ simulation Model B. Blue 
indicates regions which are bound when the kinetic, thermal and 
magnetic energies are included, green, the kinetic and thermal, 
and red, just kinetic. If only the kinetic energy is used to de- 
termine whether regions are bound, the clumps are clearly more 
extended. Two large and massive complexes are moving away 
from the spiral arms, which have high density and low velocity 
dispersions, so are strongly bound. There is also a considerable 
number of bound clumps in the interarm regions. 



regions. Overall in the bound gas there is no particular evi- 
dence of spiral structure other than the two massive 'GMCs'. 
The degree of bound gas in the interarm regions reflects the 
fact that the velocity dispersions of clumps in the interarm 
regions are lower than in the spiral arms. In the interarm 
regions there are fewer collisions and fewer interactions be- 
tween clumps. 

From hereon, we shall generally assume that the cal- 
culation of bound gas includes the thermal and magnetic 
energy, as well as the kinetic. This assumption minimises 
the amount of 'bound' gas present in the simulations. Using 
the kinetic energy alone tends to produce an unrealistically 
high fraction of bound gas. In Fig. [2] we plot the cumu- 
lative surface density of bound gas versus number density 
for model B. Here the cumulative mean surface density of 
bound gas is defined as Ebound('T') = -A^bound(w)/y4, where 
Mbound(?i) is the mass of bound gas in clumps with maxi- 
mum density > n, and A is the area of the galactic disc. It 
can be seen that the maximum density in a clump typically 
must be n > 10'' cm~'^ to obtain bound regions. This con- 
dition is less stringent if thermal and/or magnetic energies 
are ignored. Fig. [2] also shows that the bound gas represents 
6 % of all the gas in the disc. 



2.2 Star formation rate 

In order to determine the star formation rate, having deter- 
mined an estimate of the amount of bound gas, we require 
an estimate of the local timescale for star formation. As 
we mentioned above we take this to be the local dynamical 




10' 10° 10° 
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Figure 2. The mass of bound gas is plotted versus number den- 
sity n when S = 8M0pc~^ (Model B). Particles are sorted by 
density, and so n represents the maximum density in a given 
bound region. The mass indicated on the j/-axis is cumulative, 
and thus represents the total mass of bound gas in regions with a 
maximum density > n. The mass of bound gas is calculated using 
the kinetic, thermal and magnetic energy. Clumps with densities 
< 10* cm~^ do not tend to contain bound gas. The right hand 
y-axis indicates the percentage of the total gas in the disc which 
is bound. 



timescale of the gas. 



^dyn 



16G{p) 



(1) 



Here (p) is the median density of a clump. An alternative 
would be to use the volume average density, p, but we find 
that this does not lead to a noticeable difference in the re- 
sults. The star formation rate of a particular clump is then 
assumed to be M — A/bound /idyn where Mbound is the mass 
of bound gas in a clump and tdyn the dynamical timescale 
of that clump given by equation [l] The star formation rate 
per unit area, averaged over a given region, is then 



(2) 



where the summation is over all (A'') bound clumps, and A 
is the area of the region under consideration. 

It is well known that star formation is a relatively inef- 
ficient process in that not all the bound gas is converted to 
stars on a local dynamical timescale. We can allow for this 
by including an efficiency parameter e. In this case the star 
formation rate would be 



(3) 



where e is the constant star formation efficiency. However 
we do not consider the evolution of the gas once it becomes 
bound. Rather than assume a value of e for the star forma- 
tion rate, we consider what value of e would be required to 
fit observations. As we expect the value of e turns out to 
be small, justifying a posteriori our decision not to remove 
mass from the ISM in order to model star formation. 
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3 RESULTS 

We now apply these simple assumptions to the numerical 
simulations. 



3.1 Global star formation rate versus mean 
surface density 

In Fig. |3j we plot the star formation rate per unit area av- 
eraged over for the whole disc in Models A-D, estimated 
by using Eqn [2] This is then the equivalent of the global 
Kennicutt-Schmidt law. Points are shown at three different 
time frames for each surface density to illustrate the point 
that the ISM is not evolving significantly through the simu- 
lations. The points do show a correlation of star formation 
rate with surface density, roughly of the form Esfr oc E^ ®. 
However, even though the range in surface densities is not 
large, the figure does indicate that the relation we find is not 
a simple power law, with the dependence becoming shallower 
at higher densities, or equivalently showing a sharp down- 
turn at low surface densities. 

The estimated star formation rates in our models are 
much higher than observed since the formula in equation [2] 
assumes a star formation efficiency of 100%. We can however 
estimate what efficiency would correspond to observations. 
Comparing the star for mation rate at a surface density of 
10 M0 pc"^ to Fig. 2 of|Kennicutt| |l998|, we find that we 

0.05) in order to 



need to assume an efficiency of 5% (e 
produce star formation rates compatible with observations. 

We can also consider the star formation rate efficiency 
by calculating the ratio of Ebound/Etotai, as shown on Fig. 2. 
This fraction is 1.5 %, 6 %, 5 % and 6 % with surface densi- 
ties of 4, 8, 16 and 20 M© pc~^. Thus the proportion of gas 
undergoing star formation is similar in the regime where the 
slope of Eatr starts to depend linearly with E, but somewhat 
lower at the lowest surface density. 



3.2 Local star formation rate versus local surface 
density 

With effectively only 4 data points, evaluating the depen- 
dence of the star formation rate per unit area on surface 
densities averaged over the entire disc is limited. We have 
therefore divided each galaxy into 500 x 500 pc squares (the 
resolution of the recent THINGS results). We then calcu- 
late the star formation rate per unit area from the surface 
density of bound gas in each square, also assuming a star 
formation efficiency of e = 0.05, in order to compare with 
observations. In Fig.[4j we have combined the star formation 
rates from a quarter of the disc in the calculations of Mod- 
els A, B and D, which have mean surface densities of 4, 8, 
and 20 M0 pc~^ into a single plot, thus acquiring more data 
points and a greater range of densities. In the figure we plot 
for each square the estimated star formation rate per unit 
area versus the mean ISM surface density (cool and warm 
gas). 

In Fig.|4]we have overplotted our star formation rates on 
Fig. 8 of ,Bigiel et al.| ( [2008; ), which shows the star formation 
rate over the spirals in their sample versus surface density, 
as well as the globally averaged star formation rate for the 
galaxies in the survey by |Kennicxitt] ( |1998[ ) . Our simplified 
model seems to fit the observational data reasonably well. 
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Figure 3. The estimated global star formation rate 
(Sbound/^dyn) is shown versus mean surface density for the 
galaxies in Models A, B, C and D. The star formation rate 
does not show a good linear fit with surface density, and at low 
densities at least, is steeper than the observed Kennicutt-Schmidt 
law. 



Our distribution appears to support the decline in the star 
formation rate below around E = IOMqPc"'^, as seen in 
IBigiel et al. ( 2008 1 . Our points also agree with the lower 



Kennicutt (19981, similarly indicating a 



end of those from 
possible relative decline in star formation at low densities. 
In fact the star formation rate appears to show a change in 
slope, and there is much more scatter at surface densities less 
than 10 Mq pc~^. To clarify this we have binned the points 
from our results in Figure [5] There we show the mean star 
formation rate at each surface density, together with scatter 
indicated as error bars at one standard deviation. There is 
evidently more scatter at low densities and a steeper slope. 
We discuss this more in the next sections. 

Fig. [5] also shows the difference when we take only the 
regions with bound gas (blue points), and when all 500 x 
500 pc regions covering the disc are considered, even if they 
contain no bound gas (red points). For the latter we assume 
a maximum star formation rate of 10~® M0kpc~^yr^^ in 
the squares where there is no bound gas, rather than zero, 
in order to calculate the error bars. This estimate is approx- 
imately the lowest resolvable star formation rate in our sim- 
ulations. As expected the slope is steeper when regions with 
this minimal star formation rate are included (red points). 
Also the figure indicates that regions which do not contain 
bound gas tend to have average surface densities of less than 
~ 10 Mq pc-2. 



3.3 The local Kennicutt-Schmidt law for different 
tracers 

It has become apparent from recent observations that the 
star formation law varies for different tracers. The star for- 
mation law is shallower for high density gas, e.g. H2, HCN 
( |Gao fc Solomon||2004b| [Bigiel et al^|2008[ ), so Esfr goes 



roug hlyasE(l/2)' ". By comparison the star formation rate 
varies as approximately E^ ' when E is the surface density 
of all gas, i.e. Hi and H2. The law becomes even steeper 
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Figure 4. The star formation rate is plotted for the 4 (dark green 
crosses, Model A), 8 (blue crosses, Model B) and 2OM0pc~-^ 
(cyan crosses, Model D) models, assuming a star formation ef- 
ficiency of e = 0.05. The black circles are normal spirals, and 
the black squares starburst galaxies from the survey of |Kennicutt| 
( |1989|. These are all overplotted on Fig. 8 (middle right panel) 
of |Bigiel et al.| | |2008| l, which is the green and orange contour 
plot visible in the background. This shows the star formation 
rate versus surface density of Hi and H2 sampled over 500 pc^ re- 
gions from 7 spiral galaxies. The green, orange, red and magenta 
contours indicate 1, 2, 5 and 10 data points respectively. The 
star formation rates from the simulations and the observed data 
agree reasonably, and both show a downturn in slope at around 
S = 10 M0 pc-2. 



when just Hi is considered (iHcycr ct al."2004 Kennicutt & 



et. al.||2007| [Bigiel et a"Ll[2008 ), sugge sting that Hi is not a 



good measure of local star formation. |Krumholz fc Thomp-| 
(20071 propose that the transition from superlinear to 



linear occurs when the density at which a molecule is excited 
is similar to the median density of the galaxy. 

A common feature in the results from the simulations 
and the results of [Bigiel et al.|(2008| l is that there is much 
more scatter at low surface densities compared to high. A 
likely explanation is that at low densities the gas can ex- 
hibit a range of distributions - the gas can lie in a few dense 
bound clumps undergoing star formation, or alternatively 
in a more diffuse medium with very little star formation. At 
higher densities, much more of the gas in a given region is 
likely to be in bound clumps undergoing star formation and 
the possibility of larger volumes of diffuse gas is diminished. 



Bigiel et al. ( 2008 1 express this as the filling factor of the gas. 



i.e. the ratio of gas as high densities where stars are forming. 
The high density tracers preferentially select the top (high- 
est star formation rate) points from the distribution for all 
the gas. 

We can test directly the likely effect of using different 
tracers in the simulations. As we have mentioned above the 
cool component of our model ISM can be seen as a proxy for 
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Figure 5. The star formation rate for points in the simulations 
with 4, 8 and 20 Mq pc~^ (as shown in Fig. 4) is binned according 
to surface density. The scatter in the average star formation rate 
is then shown as 1-cr error bars. The different points indicate 
whether regions without bound gas are included. The blue data 
points exclude regions in which none of the gas is bound. The 
red data points are inclusive of those regions which contain no 
bound gas (and are therefore not seen on Fig. 4), where instead 
we assume a maximum star formation rate of 10^® M0 kpc^'^ 
yr-i. 



molecular gas, whereas the warm component can be seen as a 
proxy for atomic gas. It is the cool gas which is more severely 
affected by the spiral structure, and therefore the cool gas 
which is more likely to be contained in bound entities. The 
warm gas, with sound speeds comparable to the potential 
depths of the spiral potential, tends to be less affected by the 
spiral arms and less easily assimilated into bound clumps. 
Thus, although the evolution of different chemical species is 
not followed in these simulations, we can use a density cut 
to select particles over a given density. In the observations, 
it is not whether or not the gas is molecular that defines the 
slope of the power law, it is merely that the molecular gas 
traces the denser parts of the ISM. 

Fig. [6] shows the star formation rate per unit area 
against surface density for Model B when the surface density 
of all the gas is used, compared to when there is a density 
cut of 10~^^ g cm""^. With the density cut, the surface den- 
sities are only calculated using the gas with density above 
this threshold, whilst the star formation rate is the same, 
thus points are shifted to the left in the Figure. At low sur- 
face densities, much less of the gas is at high densities, so 
the points are shifted much further. Hence the slope is shal- 
lower compared to the relation for all the gas and indeed 
tends towards the observed linear relation between star for- 
mation rate and surface density for the densest gas. We also 
selected gas below this threshold, and as expected a steeper 



relation ensues. Gnedin et al. (20081 obtain similar results 



by plotting the star formation rate against densities of Hi, 
H2 and H1-I-H2 calculated in their simulations. 

In our models the slope of Espr changes continuously 
from a very steep slope at low density criterion to linear 
with a higher density criteria. The roughly linear relation 
to surface density for gas above 10~^^gcm~^ suggests gas 
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suggested by Krumholz & Thompson (20071, might be that 
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Figure 6. The star formation rate is plotted against surface den- 
sity for 500 pc^ regions from Model B, wliere tlie mean surface 
density is 8Mq pc~-^. The red crosses represent points where no 
density cut is imposed. The blue triangles correspond to points 
which only include gas of density above 10~^^ gcm~^. The slope 
is clearly shallower when only the denser material is used. Finally 



we selected gas with densities < 10 
which the slope becomes steeper. 



'gem ^ (green circles), for 



above this density is involved in star formation. Actually, 
only about a quarter of this gas is gravitationally bound. We 
therefore also calculated Esfr with a cut of 10~^^gcm~^, 
in which case nearly all the gas is bound. The slope is ap- 
proximately 1.0 in both cases, but the points are shifted to 
lower surface densities with the higher surface density cut. 

3.4 Theoretical interpretation 

There seems to be evidence from the observations, that 
for molecular tracers which presumably correspond to the 
higher density gas, Esfr oc E" with n — 1. This relation can 
be reproduced in our simulations by considering only the 
highest density gas. This relation also holds at the higher 
surface densities even when E is the total surface density. 
This appears to be approximately when E > IOMqpc"^, 
where most of the gas is cool/molecular ( jPobbs et al.|2008{ 
Krumholz et al.||2008[ ). 

This raises two questions: why is the relationship lin- 
ear, and why is it much steeper than linear at lower surface 
densities? 



3.4-1 A shallower local Kennicutt-Schmidt relation? 

At high surface densities a shallower relation of the form 
Esfr oc E^'" is found. This contradicts the straightfor- 
ward expectation of most theoretical models (see Section[L2j 
and, for example, the discussion in Section 5 of Kennicutt, 
1998), in which the star formation rate is presumed to de- 
pend on the local surface density divided by an appropriate 
timescale. For example one might take Esfr oc E/tdyn- Then 
assuming that tdyn oc E""'^, one finds that Esfr oc E^'^, in 
rough agreement with the Kennicutt-Schmidt law. 

One possible explanation for this discrepancy, as also 



star formation only takes place in the densest gas, regardless 
of the average surface density. Thus for determining the star 
formation rate, p is effectively constant, and therefore the 
local dynamical time given by a formula such as Equation [l] 
is also roughly constant. 

To examine this further, we plot in Figure[7]the various 
dynamical times of the various bound clumps against the 
mean local surface density of the 500 x 500 pc^ areas in 
which they are found in Model B, the 8 M© pc"'^ simulation. 
As can be seen from the figure, the dynamical time shows 
no particular correlation with surface density. Thus to a first 
approximation, the local star formation rate per unit area 
just depends in a linear fashion on the local mean surface 
density of bound gas. The exception in our calculations is 
the lowest surface density case (model A) where there is a 
large increase in the dynamical time at low surface densities. 
Our hypothesis differs mainly from [Krumholz fc| 



Thompson (120071) in that we do not suppose a E'"^ rela- 



tion for the total gas. Instead, as expressed in Section 3.3, 
we expect a large degree of scatter between the star forma- 
tion rate and E at low surface densities, depending on the 
local environment of the gas, whilst a linear relation prevails 
at high surface densities. 

3.4.2 T, (total) versus T,(H2) 

The difference between the tracers can be illustrated in the 
following manner. Supposing that the fraction of the ISM 
which is in molecular form, f{H2), is a monotonically in- 
creasing function of the local surface density E. As an ex- 
ample we take 



fiH2 



(4) 



for E < Eo and for some a > 0, and / = 1 otherwise. When 
E = Eo let the star formation rate be Esfr = Esfro- 

We then assume the observed relation for H2, i.e. 
Esfr oc E(7J"2) to hold at all times, and use this to cal- 
culate the equivalent relation for the total surface density 
(II1-I-II2) gas and for the Hi surface density. Clearly then for 
E > Eo, we have the linear relation Esfr oc E. 

However for E < Eo our assumptions imply that the 
star formation law for the total surface density (Hi and H2) 
is Esfr oc f{H2)T. oc T.^^'+^K Thus the relation for the total 
gas is expected to be steeper at low densities compared to 
high. 

Now consider the relation of star formation rate with 
surface density in Hi. We then expect there no longer to 
be a one-to-one correspondence between star formation rate 
and surface density! For example we now have that 



EiHI) = E-E(//2). 



(5) 



Then in our simple model (equation |4| when T,{HI) — 0, 
either there is no gas at all so that E = and Esfr = or 
the gas is fully molecular so that E Eq and Esfr J? Esfro- 
Thus overall the star formation rate is multivalued for 
a given density of Hi. The Hi surface density reaches a max- 
imum at a particular value of E, when the gas starts to 
become predominantly molecular. Above this, E(i/J) de- 
creases, but the star formation rate continues to increase. 
This is essentially the top of a parabola-like curve in T,{HI) 
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Figure 7. The dynamical time of individual bound entities is 
plotted against the surface density for the 500 pc^ regions in 
which they lie from Model B, where the mean surface density 
is 8Mq pc~^. Thus although the dynamical time of each bound 
clump is related to the density of a localised clump of bound gas, 
this relation largely disappears when considering the dynamical 
times of all the bound gas clumps in a 500 pc^ region. 



versus Espr space, and can be seen in the Hi plots of |Bigiel] 
et al.|(|2008|) and iKennicutt & et. al.l (|2007|). 



3.4.3 A change in slope 

For plots of the total surf ace density, we fin d a cha nge in the 



Krumholz et al. 



(20081 interpret 



slope around 10 M© pc 
this as the surface density at which gas becomes molecular. 
In our models, this surface density corresponds roughly to 
the density at which the surface densities of bound and un- 
bound gas are approximately equal. Below this density, the 
surface density is dominated by low density, unbound gas, 
and follows a steeper slope, as indicated by the low density 
criterion on Fig. 6. Above this density, the gas is predom- 
inantly bound, thus follows the shallower path. Since the 
gas generally is self-gravitating at densities high enough for 
molecular gas to be observable ( Hartmann et al!]|2001 l, our 
critical density is similar to that of [Krumholz et al.| ( |2008[ ). 

3.4.4 Higher surface densities 

A major disadvantage of the analysis presented in this pa- 
per is that, for computational reasons, we are unable to 
include very high surface densities, with no points over 
100 Mq pc~^. Thus, for example, unlike previous theo- 



ries and observations ( Kennicutt 1989 , Krumholz & McKee 
|2005[ ), we are not able to claim applicability of our findings 
to starbursts. 

Our simple analysis in Section 3.2.4 suggests that the 
star formation rate depends linearly on local surface den- 
sity, once the gas becomes fully molecular. Whilst a linear 
relation is observed for molecular gas measured using CO in 
normal galaxies ( [Bigiel et al.|200^ , it is not clear that this 
is the case for starbursts. Although Gao & Solomon ( 2004b I 



find a linear relation between star formation rate and lumi- 
nosity in starbursts for L{HCN), they observe that the star 
formation rate varies with CO flux as //(CO)^ "*. 

Nevertheless, we speculate here that our simple ideas 
might even be relevant to the high gas surface densities, 
provided that for these HCN is a better indicator of molec- 



ular gas mass than CO. Indeed, various authors (e.g Gao 
fc Solomon||2004a[ |Wu et al.|[2005| |Aalto| |2008^) have noted 



that CO is not a particularly good tracer of star formation 
in starbursts compared to HCN. We give two reasons why 
this might be so. 

First, the ^^CO predominantly traces warm molecular 
envelopes surrounding cold clouds ( Meier et al.|2000 Glenn 



& Hunter 2001 Gao & Solomon 2004b I. The ISM in star- 



bursts is also considerably more turbulent than that in nor- 
mal galaxies (Aalto et al. 11 1995 1. Thus it is not implausible 



to suppose that a substantially smaller fraction of the CO 
represents self-gravitating gas. Compared to our interpreta- 



tion of the Bigiel et al. ( 2008 1 results in Section 3.4.2, it may 



be that for starbursts we should regard CO as a proxy for 
HI and HCN as a proxy for H2 . 

Second, it seems plausible to assume expect that at high 
surface densities optical depth effects start to undermine 
a simple relation between L(CO) and L(HCN) and surface 
density. However, it is the CO observations which are likely 
to be effected first. Thus it may be that at high surface den- 
sities, L(CO) underestimates the mass relative to L(HCN). 

3.5 Star formation rate versus shock strength 

If, as suggested by |Roberts| ( |1969[ ) , star formation is trig- 
gered by spiral density waves, a higher degree of star for- 
mation might be expected in galaxies with stronger shocks. 
This has been the subject of much past debate. |Elmegreen| 
& Elmegreen ( 1986 1 argued that since grand design spirals 



do not show an increased star formation rate compared to 
flocculent galaxies, spiral triggering of star formation is not 
significant. Instead the gas and thus the star formation is 
merely arranged into spiral arms. Though some recent ob- 
servations show that there is a correlation between the SFR 



and spiral shock strength Seigar & James ( 2002 1 



For a stronger spiral potential, the stronger shock leads 
to more gas in the spiral arms and higher gas densities in 
the spiral shock. Therefore we may expect that the amount 
of bound gas increases with shock strength. Fig. [S] shows 
the star formation rate versus the strength of the spiral po- 
tential. The star formation rate does not show an increase 
with shock strength (within a factor of 2 or 3), but instead 
remains fairly flat. Thus a stronger shock does not appear 
to produce a higher star formation rate (averaged over the 
disc). 

The amount of bound gas depends primarily on the den- 
sity of the gas and on the velocity dispersion. Thus a possible 
explanation for the apparently constant star formation rate 
is that the kinetic energy of the dense gas also increases with 
the strength of the shock. To test this, we plot the velocity 
dispersion of the clumps in Fig. [9] against the mass of the 
clumps. The clumps in the higher shock models clearly have 
a higher velocity dispersion. Thus although they have higher 
densities, there is not a substantially greater mass of bound 
gas. However the bound gas is more concentrated to the spi- 
ral arms at higher shock strengths, and therefore there is a 
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Figure 8. The estimated star formation rate is shown versus the 
strength of the spiral potential when S = 4Mq pc~^ (Models A, 
E, F and G). The star formation remains fairly constant as the 
strength of the potential increases. 
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Figure 10. The velocity dispersion is plotted against azimuth 
for the calculation with different shock strengths. The velocity 
dispersion becomes greater in the spiral shock as the strength of 
the shock increases. 



100 




Figure 9. The velocity dispersion is shown against mass for each 
bound clump, from the simulations with different shock strengths. 
The velocity dispersions exhibit higher values at higher shock 
strengths. Best fit lines are also shown through each set of points 
to illustrate that they are offset. 



somewhat higher star formation rate in the spiral arms at 



higher shock strengths (see Section 3.6 1. 

We illustrate the velocity dispersion increase further in 
Fig. |10| where the velocity dispersion is plotted against az- 
imuth. The particles used to calculate the dispersion are 
selected from a ring of width 200 pc at a radius of 7.5 kpc. 
The velocity dispersion of the gas in the spiral arms gener- 
ally increases as the shock becomes stronger. 



3.6 Star formation in spiral arm and interarm 
regions 

As described above, the star formation rate, or mass of 
bound gas, increases with surface density, but does not vary 
significantly with spiral shock strength. Here we investigate 
whether the degree of star formation in the spiral arms com- 
pared to inter-arm regions varies according to surface den- 



sity or shock strength. Gas within a 1 kpc wide extent cov- 
ering the spiral arms is assumed to be spiral arm material. 
Figure [TT] shows the percentage of bound gas in the spiral 
arms versus surface density (blue crosses) and the strength 
of the potential (red diamonds). In these simulations, be- 
tween 65 and 90 per cent of the bound gas is located in the 
spiral arms. The percentage of bound gas in the spiral arms 
decreases with surface density. This indicates that at lower 
mean surface densities, the self-gravity of the gas becomes 
more important compared to the strength of the shock. Pos- 
sibly at high surface densities, gravitational instabilities lead 
to bound gas in the interarm regions. However the greatest 
contribution to the interarm bound gas is from gas which 
has become bound in the spiral arms, and remains mainly 
bound in the interarm region. 

The percentage of bound gas which is located in the 
arms increases with spiral shock strength. Thus a higher 
degree of star formation is likely to occur in the spiral arms 
for the models with a stronger spiral potential. Though as 
discussed in the previous section, the total star formation 
rate does not change significantly, instead as suggested by 
[Elmegreen fc Elmegreen |(1986), there is more star formation 
in the spiral arms simply because more of the gas is there. 

The mass of bound gas is calculated using the kinetic, 
thermal and magnetic energies. If only the kinetic energy 
is used, comparatively more bound gas lies in the interarm 
regions. However the trends shown on Fig.[Tl]do not change. 



4 CONCLUSIONS 

We have investigated the degree to which the relationship 
between the star formation rate per unit area and the surface 
density of the ISM in a star-forming galaxy can be under- 
stood in terms of simple input assumptions. 

We model the ISM as a self-gravitating, two-phase 
medium, with one half the mass at a fixed cool tempera- 
ture of T = lOOK, as a proxy for H2, and the other half 
at a fixed warm temperature T = lO'^K, as a proxy for Hi. 



We estimate (Section 2.2 1 the local star formation rate in 
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Figure 11. The percentage of bound gas which lies in the spiral 
arms is plotted against surface density (blue crosses) and against 
strength of the potential (red diamonds). At higher surface densi- 
ties, self gravity becomes more important relative to the strength 
of the shock, and less of the bound gas lies in the spiral arms. 
As the strength of the potential increases though (and thus the 
strength of the spiral shock), more of the bound gas lies in the 
spiral arms. 



a bound clump as being the mass of the clump divided by 
its dynamical timescale, multiplied by an efficiency factor e 
which we take to be e = 0.05 to give a fit to the observations. 

Using these simple input assumptions we find that we 
can reproduce the observed relationship which indicates that 
the local star formation rate per unit area is linearly propor- 
tional to the local surface density of dense gas (Figure |6| . 
We show that this direct linear proportionality comes about 
because (Figure [7| the local dynamical timescales of bound 
entities do not correlate with local mean surface densities. 

We also show, in agreement with the observations, that 
the total surface density (being a proxy for T,{H I + H2)) and 
the surface density of warm gas (being a proxy for E (///)) 
are not good indicators of local star formation rates. A sim- 



ple model (Section 3.4.2 1 for the dependence of H2 fraction, 
f{H2), as a function of total surface density, in the regime 
where /(-H'2) < 1, provides a simple explanation of why the 
surface star formation rate is a steeper function of total sur- 
face density in this regime. Moreover, it is evident from this 
simple model that the star formation rate does not have 
a one-to-one relationship with Hi surface density, implying 
that any attempt to correlate star formation rates with Hi 
surface density is likely to result in large scatter. 

We are also able to demonstrate from these simple con- 
siderations that the star formation rate averaged over the 
galaxy disc does not depend significantly on the strength 
of the spiral shock (Figure |8|. This is because, although a 
stronger shock does result in a higher gas density down- 
stream of the shock, it also results in a higher dispersion 
velocity (Figures |9] and 10 1. A stronger shock does, how- 
ever, result in more of the bound gas, and therefore more of 
the star formation lying close to the spiral arms. 



ACKNOWLEDGMENTS 

CLD thanks Frank Bigiel for providing an early draft of his 
paper, and the plot for Fig. 4. We also thank an anony- 
mous referee for suggestions which improved the paper. The 
calculations reported here were performed using the Uni- 
versity of Exeter's SGI Altix ICE 8200 supercomputer. This 
work, conducted as part of the award 'The formation of stars 
and planets: Radiation hydrodynamical and magnetohydro- 
dynamical simulations' made under the European Heads of 
Research Councils and European Science Foundation EU- 
RYI (European Young Investigator) Awards scheme, was 
supported by funds from the Participating Organisations of 
EURYI and the EC Sixth Framework Programme. 



REFERENCES 



Johansson L. E. B. 



Mayer L., 



Aalto S., 2008, Ap&SS, 313, 273 

Aalto S., Booth R. S., Black J. H., 
1995, A&A, 300, 369 

Agertz O., Lake C, Teyssier R., Moore B. 
Romeo A. B., 2009, MNRAS, 392, 294 

Bate M., 1995, PhD thesis, Univ. Cambridge 

Bigiel F., Leroy A., Waher F., Brinks E., de Blok W. J. G., 
Madore B., Thornley M. D., 2008, AJ, 136, 2846 

Binney J., Tremaine S., 1987, Galactic dynamics. Prince- 
ton, NJ, Princeton University Press, 1987, 747 p. 

Booth C. M., Theuns T., Okamoto T., 2007, MNRAS, 376, 
1588 

Bottema R., 2003, MNRAS, 344, 358 

CaldweU J. A. R., Ostriker J. P., 1981, ApJ, 251, 61 

Cox D. P., Gomez G. C, 2002, ApJS, 142, 261 

Dobbs C. L., 2008, MNRAS, 391, 844 

Dobbs C. L., Glover S. C. O., Clark P. C, Klessen R. S., 

2008, MNRAS, 389, 1097 
Elmegreen B. G., 2002, ApJ, 577, 206 
Elmegreen B. G., Elmegreen D. M., 1986, ApJ, 311, 554 
Gao Y., Solomon P. M., 2004a, ApJS, 152, 63 
Gao Y., Solomon P. M., 2004b, ApJ, 606, 271 
Glenn J., Hunter T. R., 2001, ApJS, 135, 177 
Gnedin N. Y., Tassis K., Kravtsov A. V., 2008, ArXiv e- 

prints 

Hartmann L., Ballesteros-Paredes J., Bergin E. A., 2001, 
ApJ, 562, 852 

Heyer M. H., CorbeUi E., Schneider S. E., Young J. S., 

2004, ApJ, 602, 723 
Hon J. L., Prantzos N., Boissier S., 2000, A&A, 362, 921 
Kennicutt R. C, 1989, ApJ, 344, 685 
Kennicutt R. C, 1998, ARA&A, 36, 189 

C, et. al. 2007, ApJ, 671, 333 
McKee C. F., 2005, ApJ, 630, 250 
McKee C. F., Tumhnson J., 2008, ApJ, 



Kennicutt Jr. R. 
Krumholz M. R. 
Krumholz M. R, 

689, 865 
Krumholz M. R. 
Krumholz M. R. 



Tan J. C, 2007, ApJ, 654, 304 
Thompson T. A., 2007, ApJ, 669, 289 
Li Y., Mac Low M.-M., Klessen R. S., 2006, ApJ, 639, 879 
Meier D. S., Turner J. L., Hurt R. L., 2000, ApJ, 531, 200 
Roberts W. W., 1969, ApJ, 158, 123 
Robertson B. E., Kravtsov A. V., 2008, ApJ, 680, 1083 
Saitoh T. R., Daisaka H., Kokubo E., Makino J., Okamoto 
T., Tomisaka K., Wada K., Yoshida N., 2008, PAS J, 60, 
667 



The star formation rate and surface density 



Schmidt M., 1959, ApJ, 129, 243 

Schuster K. F., Kramer C, Hitschfeld M., Garcia-Burillo 

S., Mookerjea B., 2007, A&A, 461, 143 
Seigar M. S., James P. A., 2002, MNRAS, 337, 1113 
Shetty R., Ostrikcr E. C, 2008, ApJ, 684, 978 
Silk J., 1997, ApJ, 481, 703 
Springel V., 2000, MNRAS, 312, 859 
Springel V., Hernquist L., 2003, MNRAS, 339, 289 
Susa H., 2008, ApJ, 684, 226 
Tan J. C, 2000, ApJ, 536, 173 
Tan J. C, Silk J., Balland C, 1999, ApJ, 522, 579 
Tasker E. J., Bryan G. L., 2006, ApJ, 641, 878 
Taskcr E. J., Bryan G. L., 2008, ApJ, 673, 810 
Tasker E. J., Tan J. C, 2008, ArXiv e-prints 
van den Bosch F. C., 2000, ApJ, 530, 177 
Wong T., Blitz L., 2002, ApJ, 569, 157 
Wu J., Evans II N. J., Gao Y., Solomon P. M., Shirley 

Y. L., Vanden Bout P. A., 2005, ApJL, 635, L173 
Wyse R. F. G, 1986, ApJL, 311, L41 
Wyse R. F. G., Silk J., 1989, ApJ, 339, 700 



This paper has been typeset from a MJi^ file prepared 
by the author. 



